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Abstract 

In this note we examine whether spherically symmetric solutions 
in Covariant Horava Lifshitz Gravity can reproduce Newton's Law in 
the IR limit A — >• 1. We adopt the position that the auxiliary field 
A is independent of the space-time metric [lOl E], and we assume, 
as in [5], that A is a running coupling constant. We show that under 
these assumptions, spherically symmetric solutions fail to restore the 
standard Newtonian physics in the IR limit A — > 1, unless A does 
not run, and has the fixed value A = 1. Finally, we comment on the 
Horava and Melby Thompson approach [5] in which A is assumed as 
a part of the space-time metric in the IR. 
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1 Introduction 



A recent power-counting renormalizable model for Gravity was proposed by 
Horava [lj. This scenario is based on an anisotropy between space and time 
coordinates, which is expressed via the scalings t — > b z t and x — > bx, where 
z is a dynamical critical exponent. Same recent papers on this topic can be 
found in Ref. [2]. 

Although Horava Lifshitz (HL) Gravity violates local Lorentz invariance 
in the UV, General Relativity (GR) is expected to be recovered in the in- 
frared IR limit. This implies a very special renormalization group flow for 
the couplings of the model, in particular it is expected that the coupling A in 
the extrinsic curvature term of the action has the behavior A — > 1, i.e. that 
it flows towards its GR value. But there is no theoretical study supporting 
this specific behavior. In addition, there are several other potential incon- 
sistencies in HL Gravity which have been discussed (see for example [3] and 
references therein). More specifically, the breaking of 4D diffeomorphism 
invariance introduces an additional scalar mode which may lead to strong 
coupling problems or instabilities, and in this way prevents HL Gravity from 
fully reproducing GR in the IR limit. 

In Ref. [I] a new Covariant HL Gravity is formulated by Horava and 
Melby-Thompson (HM), which includes an additional U(l) symmetry and 
two additional auxiliary fields A (Gauge field) and v (prepotential), aim- 
ing at resolving the above mentioned inconsistences of standard HL Gravity. 
Indeed U(l) symmetry eliminates the extra scalar mode curing the strong 
coupling problems in the IR limit, for A = 1. However, as it is shown by da 
Silva in Ref. [5], U(l) symmetry can not force the value of the parameter 
A to be equal to 1, since an action with the U(l) symmetry and A ^ 1, can 
be formulated. Note that the scalar mode is eliminated even for A 7^ 1 as 
it shown in [5J |6]. However, because the coupling A deviates from 1, stabil- 
ity and strong coupling problems (in the matter sector) arising, for details 
see [SI U\- Previous work on Covariant HL Gravity, can be obtain in Refs. 
[HI El QUI EH]- Cosmology has been examined in [5], while for spherically 
symmetric solutions, for A = 1, the reader may consult Refs. [TTJ], [TTJ . Also, 
star solutions has been obtained in Ref. [TT] . 

We would like to note that Covariant HL Gravity, as formulated by HM in 
[lj, incorporates an additional assumption for the field A, according to which 
A is assumed as a part of the metric in the IR limit, via the replacement N — > 
N — A/c 2 . Although, in the present paper we discuss the HM assumption for 
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A, we mainly adopt the alternative point of view [TQJ E], unless otherwise 
stated, according to which space time metric is given by the standard ADM 
form of Eq. (JT]) below. 

In this paper we study spherically symmetric vacuum solutions in the 
framework of Covariant HL Gravity for A ^ 1, for an action which includes all 
possible terms allowed by renormalizability requirement [12J . In particular, 
we adopt that A is a running coupling constant [5] , and we study if in the 
IR limit A — > 1 we can recover Newton's Law, which is necessary in order 
to agree with experimental data. The main result of this work, which is 
discussed also in conclusions, is that if A is a running coupling constant, 
Newton's law can not be recovered in the IR limit A — > 1. At this point it is 
worth to note that solutions for A = 1, as they presented in Ref. [10], gives 
Schwarzschild geometry expressed in Painleve-Gullstrand coordinates (for 
zero Gauge field A), so that they are compatible with Newton 's Law (and 
more generally with solar system tests). However, these solutions cannot be 
recovered continuously in the IR limit when A — > 1 (see Sec. 13.21 below). 

The paper is organized as follows. In Secf2]we summarize the most im- 
portant features of the Covariant HL Gravity. We consider then in Sec J3] the 
most general ansatz for spherically symmetric solutions, including a nonzero 
radial shift function, and we derive the equations of motion and the cor- 
responding constraints. In Sees. 13. 1[ 13.21 and 13.31 we present the solutions 
considering three situations (note that in case 13.31 only the asymptotic be- 
havior of the solutions is examined), in Sec. H]we discuss the HM point of 
view for A ,and finally Sec. [5] contains our conclusions. 

2 Covariant Horava-Lifshitz Gravity 

The action of Covariant HL gravity is structured by a set of five fields: N(t), 
Ni(x,t), gij(x,t), A(x,t) and u(x,t) (i = 1,2,3). Note that N(t), Ni(x,t), 
gij(x, t) are the standard fields that appear in the Arnowitt, Deser and Misner 
(ADM) form of the space-time metric 

ds 2 = -c 2 N 2 dt 2 + gij (dx i + N*dt) (dx j + N j dt) , (1) 

where c is the velocity of light, with dimension [c] — z — 1. In addition, N and 
Ni are the "lapse" and "shift" functions which are used in general relativity in 
order to split space-time dimensions, and is the spatial metric of signature 
(+,+,+). Note that here we are interested for the projectable version of the 
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model which implies that the lapse function N(t) depends only on the time 
coordinate. For the dimensions of "lapse" and "shift" functions we obtain 
[N] = 0, [Ni] = z — 1 . The auxiliary fields A(x,t) (potential) and u(x,t) 
(prepotential) are nondynamical fields with dimensions [A] = 2z — 2, [u] = 
z — 2 . The full action of the model is formulated as 

S = J dtd d Xy/g^N [K, tj K ij -\K 2 -V + vQ ij {2K ij + V^z/)] 

+ 2(1 - X)KWv + (1 - A) (Vis) 2 —A(R — 2SI) },(2) 

in which d = 3 is the spatial dimension, k 2 is an overall coupling constant 
with dimension [k 2 ] = z — d, and the extrinsic curvature is 

Kij ^-{<lr, V,,V, • V r V,}. i,j = 1,2,3, (3) 
where the symbol is defined as 

Qij = R ij _ }_ Rg ij + Q g ij _ ^) 

Note that this choice for z = 3 is an immediate consequence of power counting 
renormalizability request. 

For the construction of the potential term V we have considered the most 
general form which includes all possible renormalizable operators, relevant 
and marginal with dimension up to six. As the form of the potential is quite 
extended for this short report we will not present it here. For details, the 
interested reader may consult Refs. [U [10] and references therein. 

This new model, as the original HL gravity, is invariant under the foli- 
ation preserving diffeomorphism, Di f f (A4 , J 7 ) , where Ai is the spacetime 
manifold, provided with a preferred foliation structure J 7 . However, the ac- 
tion of Eq. (j2J) has an additional symmetry, in particular it remains invariant 
under a £7(1) Gauge symmetry, according to which the fields of the model 
transform as 

5 a N = 0, 5 agij = 0, 5 a Ni{x, t) = NVia 

5 a A(x, t) — a — iV'Vjai, b a v = a, (5) 

where a is an arbitrary spacetime function. Accordingly, the full symmetry of 
the action of Eq. (TSJ) is the extended Gauge symmetry: U(l) x Dif f(A4, J 7 ). 
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3 Spherically Symmetric solutions with A ^ 1 

We consider the most general static spherically symmetric metric, of the 
form: 

ds 2 = -c 2 N 2 dt 2 + -J— (dr + n(r)dtf + r 2 (d6 2 + sin 2 ddcj) 2 ) , (6) 
f{r) 

where n(r) = N r (r) is the radial component of shift functions, and N r = 
n(r)/f(r) since g rr = l//(r). 
From Eq. [2] we find the action 

S oc J dtd 3 x^^N (T — V) — Ar}, T = KijK ij - \K 2 , (7) 

in which we consider the Gauge fixing v = 0, we set Q = and A = , and 
where 



N 2 J 1 N 2 r 

r 2 .„ 2 



^i 2 / 2 + 2rW/ 2 + r/f n 2 



„ ^ = -TTf^to ('/' + /-!)■ (8) 

v/( r ) v/( r ) 

For the exact form of the potential term, see our previous work [10J. 

From variation with respect to A we find R = 0, or equivalently we obtain 

f(r) = i - 2 A (9) 

where B is an integration constant. Variation with respect to n and / give 
respectively the following two equations of motion: 

(A - 1) [r 2 nf" + 3r 2 n'f + 2 (rV + 2rn' + 2n) f] + 2(2A - l)rfn = 0,(10) 
and 

A' + y r (l-j)+^(nV-r)' = OV + (11) 



+ (1-A) 



r 4/ 
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where a prime denotes a derivative with respect to r and the differential 
operator O is defined as 



Q _rN yffN ^ n d n ( r 2 d 



%m) ■ (12) 



4f 2r 

^ n=0 

The exact form of OV can be found in our previous work [10]. Finally, the 
variation of the action with respect to N(t) gives the so called Hamiltonian 
constraint 

dr^==(T + V) = 0, (13) 

where, using a time redefinition, the lapse function N(t) is set to unity 
(N(t) = 1). 



3.1 Nonzero B ^ and zero shift function n = 

In this case, for n = , the equations of motion (TTOl) and (TTTj) and the 
Hamiltonian constraint ( TT3l) are independent of A, hence they are identical 
with the equations in the case of A = 1 which has been examined in detail in 
our previous work of Ref. [10]. We will not aim to present these results again, 
but we would like to note that spherically symmetric solutions in this case 
(n = 0) are physically relevant only in the framework HM approach when 
the auxiliary field A is considered as a part of the spacetime geometry. The 
conclusion here is that solutions, with zero shift function and B ^ 0, can be 
compatible with experiment even for A ^ 1, if we adopt HM approach for A. 
On the other hand, for the case we are interested in, where A is independent 
of the metric, this class of solutions has no physical interpretation as it can 
not reproduce Newton's Law. 



3.2 Zero B = and nonzero shift function n(r) ^ 0. 

A more interesting class of solutions is the one with nonzero shift function 
n(r) 7^ 0, for which we make the choice B = 0. Note, that if B ^ the 
potential term V blows up for r — > 0, see our previous work [10], hence, in 
order to satisfy the Hamiltonian constraint we have to introduce an unphys- 
ical lower bound in space. We would like to stress that if we set / = 1 (or 
we chose B = 0) it is possible to satisfy the Hamiltonian constraint avoiding 
this unphysical lower bound. In particular, as we will see here, the only way 
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to satisfy the Hamiltonian constrain and the equations of motion is to set 
n = 0, such that the system describes just a Minkowsky metric. 

In this case, the potential V vanishes and the Hamiltonian constraint 
reads 

/ dr{(l-A) [An 2 + Arnn + (rn') 2 ] -2(rn 2 )'} =0 . (14) 
The equation (TTOT) reads, for any A ^ 1, 



n" + -ri + —n = , (15) 



and the equation (TTTj) reads 



A' + ^W^' = [8n 2 -2rW + r 2 (n , ) 2 l . (16) 



A solution of eq. (lT5l) of the form nocr" gives a = — 1/2 ± iy7 /2, such that 

1 V7 
n(r) = — j= (acos(klnr) + 6sin(A;lnr)) , /c = — , (17) 



where a, b are constants of integration. The solution (fT7|) has to satisfy the 
Hamiltonian constraint (I14p . and one first notices that all the terms in the 
integrand are of the form 

-cos 2 (fclnr) , or -sin 2 (fclnr) , or - cos(fclnr) sin(fclnr) , (18) 

iy ry /y 

such that one can first check if it is possible to cancel the integrand. For this, 
the identification of terms proportional to cos 2 (k In r) gives 

Aabk = (1 - X)(3a/2 + bk) 2 , (19) 

the identification of terms proportional to sin 2 (ATnr) gives 

-Aabk = (1 - A)(36/2 - akf , (20) 
and the identification of terms proportional to sin(ATnr) cos(fclnr) gives 

{1 + 3X)k{b 2 - a 2 ) = {1 - X)ab . (21) 

One can easily see that the only solution for these last three equations is 
a = b = 0. 
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If one wishes to cancel the whole integral appearing in eq. fll4|) . one has to 
introduce a regularization, since the following integrals do not converge 

f°° dr f°° 

/ — cos 2 (Hnr) = / du cos 2 (ku) = oo . (22) 

J0 r J-oo 



We will therefore use instead the regularization 



L 
~k 



du cos 2 (ku) = / du sin 2 (ku) 

L J-L 

/ du cos (ku) sin (ku) = 0, (23) 



where L can be factorized in the Hamiltonian constraint (IT4|) and does not 
appear in the result. This constraint gives then 

9 -a 2 + h 2 + k 2 (a 2 + b 2 ) = 0, (24) 

such that a = b = 0. 

As a consequence, the only possibility to have A ^ 1 and f(r) = 1 is n = 0, 
which corresponds to flat space-time. On the other hand solutions, for A = 1 
and B = 0, exhibit Schwarzschild geometry, expressed in Painleve-Gullstrand 
coordinates, (for zero Gauge field A) as they are presented in Ref . [10] , hence 
they are compatible with Newton's Law, and more generally with solar sys- 
tem tests. However, we see here that these solutions (Schwarzschild geometry 
in Painleve-Gullstrand coordinates) can not be recovered continuously in the 
IR limit when A — > 1, as for A ^ 1 the model has the only trivial solutions 
with flat space-time geometry. 

3.3 Nonzero B ^ and nonzero shift function 

In the most general case when / = 1 — 2B/r (B ^ 0), Eq. (fTUj) can be 
written as 

,.2,2 2B ( „ ri n A - 2\ 
n" + -n'+-n= — In" + - + -—-) , 25 
r r z r \ 2r r z A — 1 / 

As in this case we can not find an exact result, it is reasonable to look for an 
asymptotic solution of the form 

°° f B\ p 



p=0 
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where the small parameter is B/r, and A 7^ 1 is fixed. We will discuss the 
regime of validity of this expansion at the end of this subsection. For the 
first two orders we obtain 

n 'o + l n 'o + ^ = ( 27 ) 
my , 2/nA'2/nA = 2 / r± + n A 



ry J /y \^ /y J iy2 \^ ty J ty \ ^ 27" ly^ 

such that n is given by Eq. ( TT7|) . It is then easy to see that the solution 

rii = — — ^acos(A;lnr) + 6sin(A;lnr) j , k = , (28) 



satisfies the previous differential equation, if the constants a, b are given by 

l + v7)(A-2)\ _V7 b 
4(A-1) / 2 

(l + v^7)(A-2) 
4(A-1) 

where a and 6 are the constants of integration appearing in no (see Eq. (|17p ) 
To the first order in B/r we therefore have 

I f aB \ 1 / b~B , 

n = —= a H 1 — • cos (A; In r) H — H h ■ • • smffc In r) . (30) 

'r \ r I \ r \ r 




We would like to warn the reader that, for B ^ 0, we can obtain only 
the asymptotic behavior of the solutions, as the equations of motion can not 
be solved analytically. Additionally, we have not checked if solutions with 
the above asymptotic behavior of Eq. fj30|) indeed satisfy the Hamiltonian 
constraint. However, it is clear that the oscillating behavior of Eq. ( 130]) is not 
compatible with Newton's Law: in particular we have for the potential <fi(r) = 
— n 2 /2c 2 (for details on the derivation, see [10]) the following expression 

, , . a 2 + b 2 a 2 — b 2 . , , . ab . , , s ^ / 1 \ . . 

+ cos(2fclnr) + — sin(2fc In r) + — , (31) 



2r 2r r 



where we see that the oscillating terms are of the same order of magnitude 
as the Newton potential, and cannot be canceled for nonzero values of the 
constants of integration a and b. 
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Let us come back to the expansion (TStjj) . This expansion is valid as long 
as (B/r)n p+ i << n p , and we define the critical distance r c below which the 
expansion is not valid anymore, and (B/r c )n p+ i(r c ) ~ n p (r c ). If we consider 
then the first two terms of the expansion, as well as the expressions ( 129]) for 
a and b, we obtain 

1 B 1 

W c " W A^l ' (32) 

such that r c ~ B/(X — 1). Finally, for r « r c and from Eq. (jTUI) . we see 
that n(r) — > when A — > 1, while the asymptotic behavior of Eq. (TStjj) is 
valid only for r >> r c , with r c — > oo when A — > 1. 



4 Comments on the HM interpretation for A 
as a part of spacetime metric 

The IR asymptotic behavior for A, which now is assumed as a part of space- 
time metric, can be obtain by using Eq. ffTTl) : 

9 A Kx-AqB t ^ cos(2Hnr) T ^ sin(2Hnr) n f l\ , N 
= c 2 +A + — —+^2^ ^+X 3 ^ -+0 - , 33 

y^" ,y ?" \ 7" / 

in which we have kept only the leading terms of the order of 1/r, and Aq is a 
constant of integration. In the expression for Am, the constants Ki, K 2 , K 3 
are 

K x = Q(l-X)(a 2 + b 2 ) (34) 
K 2 = h(\ + Q)(b 2 - a 2 ) + 2v / 7(5A - 4)a6] 

o 

K 3 = i[v / 7(5A-4)(6 2 -a 2 )-2(A + 6)a6] . 

o 

In what follows we set Aq = — c 2 . 

In the case of HM theory the "Newtonian" potential, if we drop higher 
order terms (1/r™ with n > 2), is given by the equation 



2c 2 

-Zl _ cos(2Hnr) ,_ sin(2/clnr) / 1 
— + Z 2 i '- + Z 3 — * ^ + O — 

ry ry ry \ 



where 

25 



Z x = [Q\ - ^-j (a 2 + b 2 )-B (36) 
Z 2 = -^[(A + 4)(6 2 -a 2 ) + 2 v / 7(5A-4)a&] 

o 

Z 3 = -i[ v / 7(5A-4)(6 2 -a 2 )-2(A + 4)a6] . 

o 

where a = a/c and b = b/c. It is possible to cancel simultaneously Z 2 and 
Z 3 only for a = b = 0, since the following system of linear equations 

(A + A)x + 2 V / 7(5A - 4)y = (37) 
v / 7(5A-4)x-2(A + 4)?/ = , 

where x = b 2 — a 2 and y — ab has negative determinant, D = — 2(A + 4) 2 — 
14(5A — 4) 2 < 0. However, if we set Zi = 2M and keep suitably small the 
parameters a, b « \/M (so small that the oscillations are not observable) 
we can recover the Newtonian potential (and the Schwarzschild metric) even 
for A 7^ 1, so HM theory can not fix A to unity. In addition for A = 1 the 
oscillating terms do not vanish. 



5 Conclusions 

We have studied spherically symmetric solutions of Covariant Horava Lifshitz 
Gravity for A ^ 1, namely we have assumed that A is a running coupling as it 
presented in Ref. [5]. We found that, for zero radial shift function n(r) = 
and B ^ 0, spherically symmetric solutions are independent on the running 
coupling A, hence they are identical with solutions when A = 1 which have 
been analyzed previously in [TO]. For nonzero radial shift function n(r) ^ 
and B = 0, we show that the only solution of the model is the flat space-time 
metric if A ^ 1. In the general case: n(r) ^ and B^O, numerical work is 
necessary for analyzing the spectrum of solutions and checking if they satisfy 
the Hamiltonian constraint, which is beyond the scope of this note. However, 
in this case it is possible to obtain the asymptotic form of the metric which 
is enough in order to compare with experimental data. 

In this paper we mainly focus to the question whether spherically sym- 
metric solutions can be compatible with experiment in the IR limit A — > 1. 
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In particular we have shown that if we adopt the position that A is indepen- 
dent of the space-time metric (hence the Newtonian potential is given by Eq. 
(j3ip ). it is impossible to recover Newtonian physics in the IR limit. If B = 
the only spherically symmetric solution has flat space-time metric, so in the 
limit A 1 we can not recover continuously solutions for A = 1, which are in 
agrement with solar system tests (see Ref. [ID])- Also, for 5^0, oscillating 
behavior of solutions in the IR limit is non compatible with Newton's Law. 

According to the altenative HM approach, A is assumed as a part of the 
potential according to the equation (|35|) . In this case, as we show in Sec. HJ 
it is possible to recover Newton Law (or more general Schwarzschild space- 
time geometry) by adjusting suitably the integration constants a,b « y/M, 
independently of the value of A. As a result the limit A — > 1 is not necessary 
in order to achieve agrement with solar system tests (or equivalently even for 
A / 1 we could have agreement with experiment). 

From the above discussion we conclude that, if we adopt that the auxiliary 
field A is not a part of the metric, Newtonian potential can not be recovered 
in the IR limit A — > 1. Therefore we have to invoke a mechanics for fixing 
A to unity (for example a new symmetry) and construct a model which is 
physically relevant, as only solutions for A = 1 and B = are compatible 
with experimental data. 
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